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The correlation functions of an arbitrary number of boundary monomers in the system 
of close-packed dimers on the square lattice are computed exactly in the scaling limit. The 
equivalence of the 2rt-point correlation functions with those of a complex free fermion is 
proved, thereby reinforcing the description of the monomer-dimer model by a conformal free 
field theory with central charge c = 1 . 

1. INTRODUCTION 

The dimer model has been originally introduced to describe physical adsorption of diatomic 
molecules on crystal surfaces [l|]. The first studies of the dimer model are dated the early sixties, 
with pioneering works by Kasteleyn Q], Fisher Q], Temperley and Fisher 4\, Ferdinand [5] and 
Wu g], who have studied the number of close-packed dimer configurations on specific (mostly 
square) finite bipartite lattices. Soon after that, correlations between dimers and monomers on 
the square lattice have been examined in [7]. The effects caused by the msertion of monomers on 
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the square lattice have been reconsidered in a number of recent works 
monomer correlations have also been revisited recently, with significantly different conclusions, on 
the triangular lattice (a non-bipartite graph) in [12]. 



A close-packed dimer configuration on a square grid is an arrangement of dimers such that a 
dimer covers two adjacent sites and every site is covered by exactly one dimer. As a generalization, 
one can look at the dimer configurations when some sites, called monomers, are not allowed to be 
covered by dimers. Corresponding to these arrangements of monomers and dimers, one considers 
the partition function 

Z(w h ,w v \zi,. . . ,z N ) = wl h Wv v - (1.1) 

coverings 

It counts the number of dimer coverings in presence of N monomers located at positions z\, . . . , zjy, 
in the bulk or on boundaries, with weights Wh and w v assigned to horizontal and vertical dimers. 
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As the number + n v of dimers is fixed, the partition function essentially depends on Wh,w v 
through the ratio Wh/w v only. 

The dimer model belongs to the class of so-called free-fermion models [lJ. It is well known 
that the partition function of the free-fermion models admits a representation in terms of fermionic 
Gaussian integrals which leads to determinant expressions for the partition and correlation func- 
tions. Despite the simple form of the lattice action in the free-fermion representation, the evaluation 
of correlation functions for certain physical observables can be more complicated because local vari- 
ables correspond often to non-local fermion correlators. As illustrated in [7(], the monomer-monomer 
correlations in the dimer model need typically non-local computations if monomers are located in 
the bulk of the lattice. However the situation simplifies significantly for boundary monomers. In 
this article, we take advantage of this simplification and compute all monomer correlations ex- 
actly, in the scaling limit. We will show that in the continuum limit, the correlation functions for 
boundary monomers can be expressed in terms of complex free fermions located at sites occupied 
by the monomers. This implies that the boundary monomer degrees of freedom are described by 
a conformal field theory with central charge c = 1. 

Our result complements previous results related to the description of the general dimer model 
by a conformal field theory. If a conformal field theory with central charge c = — 2 accounts 
well for the dimer degrees of freedom 



14(, it seems that the full monomer-dimer model should 
correspond to a conformal theory with c = 1. An early indication of this can be traced in the 
work of Au-Yang and Perk who noted a close relationship between the 2- monomer correlator 
computed by Fisher and Stephenson [7] and the squared correlator of two Ising spins, equivalently 
the correlator in a doubled Ising model. Likewise dimer correlations on the square lattice have 
been more recently reinterpreted as correlators of two uncoupled massless Majorana fermions (and 
by massive Majorana fermions in the perturbation away from the square lattice to the triangular 
lattice) by Fendley et al 12]. More generally, Kenyon [l(| has shown that a quantity associated 
with a dimer configuration, called the height function, converges in the scaling limit to a Gaussian 
free field (or free boson). Like the complex free fermion, the free boson corresponds to a conformal 
theory with central charge c = 1. The relation between the dimer model and free fermions has 
been revisited very recently by Dijkgraaf et al [l7]. 



FIG. 1: Orientation of the bonds defining the signs in the Kasteleyn matrix K. 

2. GENERAL SETTING 

Here we briefly recall the way the partition function can be calculated when the graph is a finite 
portion of a square lattice. In the following, we set Wh = w v = 1 as we are not interested in the 
directional properties of the dimer coverings. 

The simplest and best understood situation is when all sites of the domain C must covered by 
dimers, the so-called close-packed limit. In this case, the partition function, which simply counts 
the number of dimer configurations so that every site of £ is covered by one dimer, can be expressed 
as the Pfaffian of an antisymmetric matrix K, 

Zf ] = Pf K = +v / det K. (2.1) 

The sign is chosen so that the partition function is positive, and the matrix K is a weighted 
adjacency matrix of C. Each term in the expansion of the Pfaffian is naturally associated with a 
dimer covering, but the entries of K have to be suitably chosen so that each covering contributes 
1 to the partition function. 

There are different choices for K, but a convenient one for what follows is the matrix originally 
considered by Kasteleyn. K is an oriented adjacency matrix for C, so that Ky = if the sites i,j 
are not nearest neighbours, and Ky = ±1 otherwise. The signs are given pictorially in Fig.l: an 
arrow from i to j indicates that = +1 and Kji = — 1. The essential property of K is that the 
product of entries around any elementary cell of C is equal to —1, 

ii i_2 

K ili2 K i2i3 K i3i4 K i4il = -1, for any cell Q . (2.2) 

i 4 i 3 

As K is closely related to the Laplacian on £, standard techniques like Fourier series can be used 
to compute the determinant of K, and then its Pfaffian. 
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Monomers are sites which cannot be covered by dimers. Thus dimer configurations on C in 
presence of monomers at z±, Z2, ■ ■ ■ , zn are close-packed dimer coverings of C \ {zi, Z2, ■ ■ ■}■ 

If all of the N monomers are on the boundary of £, the matrix K(z\, Z2, ■ ■ •) defined from 
K by removing the rows and columns labelled by the sites Zi,z 2 , ■ ■ ■ is the Kasteleyn matrix for 
C \ {zi, Z2, ■ ■ .}. As it still satisfies the property (12, 2p . one has 

Z { c N \ Zl ,Z2, . . .) = Pf K( Zl ,Z2, ...) = + V / det K( Zl , z 2 , . . .)• (2.3) 

The monomer correlations are then defined by the ratios 

c &>**> • • •) = ^) - ■ (2 ' 4) 

We are eventually interested to compute these ratios in the thermodynamic limit. It is then much 
more convenient to express K(z%, Z2, ■ ■ •) = K + B(z\, Z2, • • ■) as a finite rank perturbation of K, 
localized around the monomer positions. If indeed B has rank iV (that is, Bij = except if i, j are 
in a set of N sites), the ratio of infinite-dimensional determinants 
Z C N (zi,z 2 ,...) 
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^det K~ l K(z u z 2 , ■■■) = V / det(I + K~ l B), (2.5) 



reduces to a finite, rank iV determinant (which however involves entries of the infinite-dimensional 
matrix K^ 1 ). 

The defect matrix B has to satisfy two requirements: (i) if k is a neighbouring site of a monomer 
located at z, then B z ^ = —K z ^ = —B^ z so that z is effectively cut off from the rest of the grid; 
and (ii) the restriction of B to the monomer sites z\ , Z2 , ■ ■ ■ must have a determinant equal to 1 (the 
simplest solution is to set B Xi Zi = 1, but this is not always the most convenient way, see Section 
4). All other entries of B are equal to 0. The rank of B increases linearly with the number of 
monomers, and so does the size of the determinant. 

When some of the monomers are away from the boundary, the situation changes dramatically. 
Removing a non-boundary site from C creates a new elementary cell, around which the product of 
the restricted K matrix elements is not equal to —1. This can be remedied by changing the signs 
of K along a path going from one monomer to another monomer. This in effect introduces a non- 
local defect matrix, and complicates the calculation since the size of the determinant increases with 
the distance between the monomers. To date, the only known exact result on the square lattice 
is the old result by Fisher and Stephenson [7], who proved that the 2-point correlation of bulk 
monomers in the scaling limit decays like r -1 / 2 . In contrast the monomer and dimer correlations 
decay exponentially on the triangular lattice |l2|. There are indications that this last behaviour 
holds on non-bipartite graphs. 
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ISOLATED MONOMERS ON A BOUNDARY 



Our purpose is to calculate the monomer correlations (12. 5p for an arbitrary number of isolated 
monomers on the boundary of the discrete upper half-plane (UHP), with no monomer away from 
the boundary. For definiteness, we take the boundary to be the line y = 1, so that the discrete UHP 
corresponds to {(x, y) : x S Z, y € Z>o}. The positions of the monomers along the boundary are 
denoted by (xi, 1), and we are interested in the scaling regime where all distances 
large. We assume the Xi are ordered from left to right, so that X{j < for i < j. 

The first ingredient we need is the inverse of K. As mentioned in the previous section, the 
matrix K itself is defined from Fig.l where the rectangle is extended to the UHP. The orientation 
of all horizontal bonds is to the right, while that of the vertical bonds alternate; we fix the reference 
point by deciding that the vertical bonds on the line x = are oriented upwards. The matrix is 
then given by 

K(x,y),{x',y') = [Sx,x'-1 ~ d x ,x'+l]8y,y' + (— l) 2 5 x ,x' [&y,y'-l ~ <W'+l]- (3-1) 

It is instructive and easy to compute the square of K, 

' number of nearest neighbours of i, for i = j, 
-1, if f-j = ±(2,0) or ±(0,2), (3.2) 



otherwise. 



One sees that — K 2 connects sites within each of the four sublattices corresponding to the parity 
of the x-coordinates and of the y-coordinates, and that its restriction to any one of these is equal 
to a specific Laplacian. On the odd-odd and even-odd sublattices, the Laplacian is subjected to a 
closed boundary condition since the sites on the boundary (y = 1) have three nearest neighbours. 
On the other two sublattices, which have the line y = 2 as boundary, along which the sites have 
four nearest neighbours, the Laplacian is subjected to the open boundary condition. We can write 

— K 2 = A£ dd)0dd © A^ cn odd © A odd even © A°P en even . (3.3) 

This allows to find the inverse of K, 



K- 1 



c cl ffi r" cl ffi c op ffi n op 

^oddjOdd ^ ^everijOdd 533 ^odd^ven ^ "evei^even 



K, (3.4) 

in terms of the well-known Green matrices G = A -1 . The Green matrices G cl and G op are related 
to the inverse Laplacian G on the full (discrete) plane by standard formulae, 

G c \x - x'; y, y') = Gf x ^ [xl , yl) = G(x -x',y- y') + G(x - x ', y + y '-l), (3.5) 
G op (x - x'- y, y') = Gf xA{x , >yl) = G(x -x',y- y') - G(x -x',y + y'). (3.6) 
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One can now write the inverse of K explicitly: 



R- 1 



G cl (^=f=i; «±i, ^±i) - G cl (2^+i; ^±1)1 , if (x - x')yy' = 1 mod 2, 
G op ( ^l . | ; _ G op ( ^+l. £ ; if (x _ _ l){y , _ 1) = 1 mod 2) 

(_!)* \ G c \^f- »±i, K^) - G cl (^; *±i, £ )1 , if (x - x' - l)y(y' - 1) = 1 mod 2, 



(-If G op ( 



i x—x' . y y'+l 
2 ' 2' 2 



) - G° P ( 



x—x' . y_ y'—l s 
2 ' 2' 2 • 



, if (x — x')(y — l)y = 1 mod 2. 



0, otherwise (when x — x' and y — y' have the same parity). 



(3.7) 



Bi(x) 



(3.8) 



By using the relations (|3.5|) and (|3.6|) . one may check that K~ l is antisymmetric. 

We are now ready to compute the relevant determinants (|2.5p . If there are iV monomers on the 
boundary, the B matrix is 4iV-dimensional. A possible choice is to write it as the direct sum of ./V 
4-by-4 blocks B\(xi), one for each monomer, so that B = ®iB\(xi), with 

-10 
{-If 

V i o o o / 

in the basis where the indices 1, 2, 3, 4 correspond respectively to the monomer itself, its left, its 
upper and its right neighbour. 

Using this explicit form of B as well as the inverse of K, one may easily compute the correlation 
functions, given by 

Z(X1,X 2 , • • • ,X N ) 



C(xi,x 2 , . . -,x N ) 



(3.9) 



For large distances \xij\ ^> 1, the asymptotic form of these correlators may be obtained by using 
the following expansions of the Green matrix on the plane, 



G(m, k) = G(m, 0) 



k 1 



k 4 - 3k 3 
Anm 2 87rm 4 



+ 



k <C m, 



G(m, 0) = log \m\ + . 

2ir 



(3.10) 
(3.11) 



We have computed explicitly the first few correlators, and found the following expressions. The 
1- and 3-point functions are identically zero, 



C(xi) = C(X1,X2,X 3 ) = 0, 



(3.12) 



as one would expect. On a finite rectangular grid (with an even number of sites), there must be an 
even number of monomers since otherwise, the rest of the rectangle cannot be covered with dimers. 
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The 2-point function is equal to, at dominant order 

2 



+ . . . , if X12 is odd, 
C{x 1 ,x 2 ) = { ™12 ' (3.13) 

0, if x\2 is even. 

Again this characteristic difference in the parity of the distance between the two monomers is 
expected, since on a finite rectangle, there must be an equal number of monomers on the even 
sublattice as on the odd sublattice. 

For alternating even and odd monomer positions 2j (that is, all Xj — Xj_j_i are odd), the 4- and 
6-point functions are equal to 

C(xi,---,x 4 ) = -U— *— + — *— ) + ••• ( 3 - 14 ) 

IT 2 I £12X34 X14X23 J 

C(xi, • • • ,2 6 ) = 

M 



r 1 1 1 1 1 , 
J + h h + } + ... 



7T 3 I £12X34256 212X36X45 X14X23256 214X25X36 Xi 6 X 2 32 4 5 2i 6 2 2 5 234 • 

(3.15) 

In the scaling limit, these correlators exactly match those of a complex chiral free fermion if), 
lim C(xi, • • • , x 2n ) = (^(21) ^(22) ^(23) ■0 t (24) • • • ip(x 2n -i) ^ t (2 2n )), (3.16) 

scaling 

if the 2-point functions are normalized as 

(V(x)^(y)) = {^{x)^{y)) = - 2 (3.17) 

vr(2 - y) 

<V>(2)V(y)> = (^(x)^(y)) =0. (3.18) 

Physically the charged fermions V = V'e and ^' = ip can be interpreted respectively as the 
insertion of a monomer at an even position and at an odd position (or vice- versa), so that a 
globally neutral correlator indicates an equal number of even and odd monomers. Their real and 
imaginary parts, ip\ = ^= {tp + ift ) and V2 = J^i^ ~ V^) can De viewed as two uncoupled Ising 
Major ana fermions which form a chiral conformal field theory with c = 1. 

We prove, in the next section, the equality (|3.16p for an arbitrary value of n. 



4. BOUNDARY MONOMER CORRELATORS 



We start with the cases where the correlators do not vanish identically: we place 2n monomers 
on the boundary, located at (xj,l), n of which are at even (odd) positions. Without loss of 
generality, one may choose xi even, X2 odd, X3 even, and so on. 
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FIG. 2: Labelling of the 8n sites involved in the calculation of the determinant. The solid circles represent the 
positions of the monomers (ari, X2, . . .), all located on the boundary, the open circles their nearest neighbours. 

The matrix B used in the previous section is not the most convenient choice to carry out the 
general calculation. We slightly modify the entries of B labelled by the monomer positions: we 
set the diagonal elements to and connect the monomers by pairs, by setting Br Xu \\i x j\ = 1 = 
— B( Xi+1 ,i),(xi,i) f° r an * odd. It means that the restriction of K + B to the monomer sites is not 
the identity matrix like in the previous section, but a direct sum of 2-by-2 blocks equal to ^ _^ 
whose determinant remains equal to 1. The other off-diagonal elements of B are as before. 

This choice ensures that B is also antisymmetric, and such that Bij = if i— j has coordinates of 
equal parities. Since the matrix K -1 has the same property, see (|3.7p . it follows that (I+K^B)^ = 
if i—j has coordinates of opposite parities. By an appropriate ordering of the site indices, I+K~ 1 B 
can thus be brought to a block-diagonal form. 

The determinant to be computed has dimension 8n since B has rank 8n: there are 2n monomer 
sites, and each of them has three nearest neighbours. We label the sites as in Fig. 2, using two 
types of roman indices, a and a, each type of label taking its values in {1, 2, • • • , 4n}. The labelling 
is such that the differences a — b or a — b have coordinates of equal parities, and that differences 
a — b or a — b of unlike sites have coordinates of opposite parities. In this basis, the matrices B and 
K^ 1 are off-diagonal, B ab = B ab = K~ b l = K~ b l = 0, while the matrix 1 + K~ l B is block-diagonal, 

((l + K^BU \ 

(I + K-'B)^ = . (4.1) 

V (I + K^BUj 

Moreover the two diagonal blocks are closely related. As is manifest in Fig. 2, the two types 
of sites are exchanged by a mirror symmetry, under which the monomer coordinates and the 
separation distances are transformed according to X£ — > X£ = x 2n +i-e and — X£ — > X£ — It 
follows that the second block depends on the separation distances xg — x^ in the same way the first 
block depends on Xk — xe, or equivalently, 

(I + K-\x k - x e )B) ab = (I + K-\x t - x k )B) ab . (4.2) 
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As it turns out, the determinant of the first block will be invariant under the substitution 
Xk — xg — > X£ — Xk, and therefore equal to the determinant of the second block. Putting all 
together, one obtains the correlations as 

C( Xl ,x 2 ,...,x 2n ) = ^/det(I + K~ l B)ij = \det(I + K- 1 B) ab \. (4.3) 

We finish the proof by showing that this last determinant reduces, in the scaling limit, to the 
fermionic 2n-point function (|3.16j) . 

We start by grouping the 4n sites a by four, each group receiving a label I, between 1 and n: 
1 = 1 corresponds to the first four sites {1, 2, 3, 4}, I = 2 to the next four sites {5, 6, 7, 8}, and so 
on. Accordingly we write the matrix (I + K~ 1 B) a b = Ajj in a block form, where all blocks Ajj 
have dimension 4. 

Because the matrix B does not connect sites belonging to different values of I, the diagonal 
blocks An are ah given in terms of a single matrix function A%, as are the off-diagonal blocks in 
terms of a second matrix A 2 . The entries of An only depend on the distance x 2 i-\ 21, so that A± 
depends on a single variable, Aji = A\ (2:21-1,2/); likewise the off-diagonal block Ajj only depends 
on the two distances a?2/-i,2J and x 2 i i2 j, so that Ajj = A 2 (x 2 i-\ >2 j;x 2 i j2 j) (to see that the other 
two distances #21-1,2,7-1 and x 2 i,2J-i do not enter in A 2 requires a simple calculation, illustrated 
below) . Therefore the matrix of which we need to compute the determinant has the form 

/ Ai{xi 2 ) A 2 {xu;x 2i ) A 2 (x 16 ;x 2e ) ■ \ 
A 2 {x?, 2 ;x i2 ) A 1 (x 34: ) A 2 (x 36 ;x i6 ) ■■ 
A 2 (x 52 ;x 62 ) A 2 (x 5 4;x e4: ) Ai(x 56 ) 

\ ; ; i ••./ 

Our next task is to compute the asymptotic value of the blocks A\ and A 2 , when the distances 
become large. This is straightforward as soon as one is familiar with the notations and with the 
form of B. Suppose that we want to compute the (l,l)-element of the off-diagonal block Ai = \^j =n 
corresponding to the two subsets of sites fully displayed in Fig. 2. The first labels in / = 1 and 
J = n are respectively the sites 1 and 4n — 3. One finds 

(4r=i,j=n)(i,i) = (I + K~ 1 B)i^ n _3 

= ^l,l^l,4n-3 + ^l,2-^2,4n-3 + -^lJ-%,4n-3 + -^l,4^4,4n-3 

= ~ K (x\,l),(x 2 „,l) - K (x\,l),(x 2n - 1 +l,l) ~ K (x\,l),(x 2n - 1 ,2) + ^(xi,l),(s2n-l-l.l)' ^' 5 ) 



(I + K~ l B) ab = Aj j 



(4.4) 
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The last three terms cancel because they are equal to (K l K )(xi,i),(x 2n -xA) = leaving, for 

X — X\ X 2n , 



(A I=ltJ=n ) n = -K(l lUx2n>1) = -G d (^; 1, 1) + G^; 1, 1) 
= -G{^, 0) - 1) + G(2±i,0) + G(2±l, 1) 



2 ' 



+ ... 



(4.6) 



for large x by using (|3.1ip . Similar calculations for the other entries and for the diagonal blocks 
yield the matrices A\ and A 2 explicitly as 



^i(x) 



/ -- 

TTX 


_2_ 

TTX 


7ra; 




1 

TT 


1-i 

7T 


_ j_ 

7T 


j_ 

7T 


1-2 

7T 


-1+2 

7T 


2 

TT 


1-2 


V -\ 


j_ 

TT 


j_ 

TT 


1-1/ 



+ 



A 2 (x;y) 



1 


2 


2 


2 


2 


\ 




7T.T 


7TX 


TTX 




__2_ 


^_ 


J_ 


__2_ 










",V 


Try 




















\ 


__2_ 


_2_ 




__2_ 


/ 


wy 


".V 







+ ... 



(4.7) 

where the dots represent lower order terms in x or y. 

A first observation is that the full matrix Ajj contains exactly n lines with all their elements 
of order —1 in the distances X£ m , while all the other 3n lines contain elements of order in 
these variables, coming from the A\ blocks. From this, it follows that the dominant term in the 
determinant has order — n in the distances. Anticipating that the coefficient of this term does not 
vanish so that the scaling dimension of the 2n-correlator is n, we may neglect the y dependence in 
the ^2-blocks, and use the simplified matrix, 



f- 


2 

TTX 


2 

TTX 


2 

TTX 


2 \ 

TTX 









































) 



A 2 {x) = A 2 {x; 00) 



At this stage, the full matrix has the following form, 

/ At(x 12 ) A 2 (xu) A 2 (x w ) 

^2(^32) ^4l(^34) -42(^36) 



(4.8) 



(l + K- l B) ab = Au 



^2(^52) ^2(^54) Ai(x m ) 



V 



(4.9) 



/ 



It has n block columns, each formed of four columns. Within each block column, we uniformly add 
the first column to the second and third ones, and subtract it from the fourth one. This column 
operation does not change the value of the determinant, and can be done at the level of the small 
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matrices A\ and A 2 . Doing this recasts the determinant into the following form, 



det(I + K~ l B) ab = det 



V 



2 


u 


n 
U 


U 


2 


n 
U 


n 
U 


n 
u 


2 


n 
U 


n 
U 


n 
U 




7TX12 


TTX14 






I 

7T 


1 


































1 - - 

7T 





1 































7T 








1 




























2 


u 


n 
U 


n 
U 


2 


n 
U 


U 


n 
U 


2 


n 
U 


n 
U 


n 
U 






34 


















7T 


1 


































1 - ^ 

7T 





1 































| 








1 
















2 











2 











2 















7ra:54 






























1 

7T 


1 


































1-2 

7T 





1 































_ 1 

IT 








1 













(4.10) 



/ 



It clearly factorizes into the product of two determinants. One has order 3n, made of all rows and 
columns except the 1st, 5th, 9th, and is equal to 1, while the other contains the remaining rows 
and columns. We therefore obtain 



lim C(xi,...,x 2n ) = hm | det(I + K~ 1 B) ab \ = ( — 

scaling scaling V TT 



2\n 



det 



1 



^2i-l — %2j' 1<M<™ 

This last form, a Cauchy determinant, can be evaluated explicitly (see for instance 

-2\n Ul<i<j<n ( x 2i-l ~ X 2 j-l){x 2 j ~ X 2i ) 



(4.11) 



W) 



lim C(xi, . . . ,x 2r , 

scaling 



TT 



(4.12) 



Ul<i,j<n ( X 2i-1 ~ X2j) 

and exactly reproduces the free fermion correlator (|3.16p . 

We complete the proof by showing that all the correlators which do not contain an equal number 
of even and odd monomers vanish identically. Without loss of generality, one may assume that in 
addition to the 2n monomers at alternatively even and odd positions (xj, 1), there are M monomers, 
all located at either even or at odd positions (y,, 1). We consider the case where the M additional 
monomers are on even sites, the other case being similar. We assume x\ < x 2 < ■ ■ ■ < x 2n < yi < 
• • • < DM- 

The matrix (I + K~ l B) has now dimension 8n + AM: the first 8n labels will be as above, 
4n indices a and 4ra indices a. The other 4M sites will be ordered in a more natural way: first 
monomer at y±, its left, upper and right neighbours, and so on for the others. The restriction of B 
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to the 8n sites is kept as above, with oriented bonds between x-monomers. As M can be odd, we 
take on the other 4M sites a direct sum of -Bi(y«) matrices, written in (j3.8j) . 

Let us now compute the n + M rows of (I + K~ 1 B) labelled by the even x-monomers (corre- 
sponding to x±,xs, . . .) and by the y-monomers (in the case where all %n are odd, one would look 
instead at the odd x- monomers). A straightforward calculation shows that the rows labelled by 
the even monomers (xi, 1), i odd, are equal to 



(I + K 1 B)( XiA) ^. - y K ( x ],l),(x 2 ,l)' K (x],l),(x 2 ,l)> K (x],l),(x 2 ,l)' K (x],l),(x 2 ,l) 5 
~ K (xi,l),(x 4 ,l)> K (xi,l),(x 4 ,l)i K (x x ,l),(x 4 ,l)> ~ K (x x ,l),(x 4 ,l) i "' '■> 



(x 



],l),(x 2 n,l)> K (x],l),(x2 n ,l)> K (x],l),(x 2n ,l)> K (x],l),(x 2n ,l) ' ''"'°; 0,---, 0^.(4.13) 



The first 4n entries are non-zero, proportional to an matrix element of K~ l ; then there is a group 
of 4n zeros corresponding to the underlined indices and another group of zeros corresponding to 
the last 4M sites. 

Because all yi are even, the first 8n entries of the rows (I + K~ 1 B) yu . are exactly given by the 
formula (|4.13p where one simply replaces X{ by yi. Moreover, one may check that the last 4M 
entries are also zero, so that one obtains the simple result that 

{l + K- 1 B) {yi ^. = {l + K- 1 B) {Xi ^. . (4.14) 

By column additions and subtractions, one may bring these rows to the form where only n 
entries are non-zero, for instance, 

(i+^-^W), = (-^;,i),(, 2 ,i)' ' ' ' -^;,i),(x 4 ,i)' o ' o '°; -^k(w)> o ' o '°; 

0, ■■■,(); 0,---,oV (4.15) 



and similarly for (1+ K~ 1 B)( y .^^.. In this form, the full matrix, which has the same determinant as 
the original one, contains n + M rows which are vectors in an M. n vector subspace. Consequently, it 
has at least M left eigenvectors with zero eigenvalue, therefore its determinant vanishes identically, 
and so does the corresponding correlator. 



5. STRING OF MONOMERS 



Up to here, we have considered isolated monomers, far apart from each other, and lying on the 
boundary of the upper-half plane. We have shown that the correlation functions for such monomer 
configurations can be understood, in the scaling limit, as free fermion correlators in a conformal 
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field theory with central charge c = 1. For completeness, we briefly discuss a different situation, and 
in some sense opposite, namely the case where the monomers form a compact cluster of consecutive 
sites. So we consider 2n monomers on consecutive boundary sites, and ask for the corresponding 
correlation Cm, defined as before as the ratio of the partition function with the 2n monomers to 
the partition function with no monomer. 

Since the 2n monomers can be covered in a unique way by n dimers, we may think of the 
dimer configurations in presence of the monomers as close-packed dimer configurations with a 
fixed string of n consecutive dimers, all oriented along the boundary. Being now formulated as 



a pure dimer problem, with a prescribed 



Temper ley correspondence (see also 
trees, can be used. 

Given a dimer configuration, an associated configuration of arrows is defined on the odd-odd 
sublattice £ dd (it contains a half of the boundary sites of the upper-half plane): if a dimer touches 
a site of £ dd) one draws an arrow from that site towards its nearest neighbour in £ dd) i n the 
direction of the dimer; those dimers which do not touch sites of £ dd are uniquely fixed once the 
dimers which do touch £ Q dd are given. The so-obtained arrow configuration has the property that 
it cannot form closed loops, because a loop would encircle an odd number of sites of the original 
lattice, which therefore could not be fully covered by dimers. So the arrow configuration defines a 
spanning tree on £ dd- 

In this correspondence, a prescribed string of n consecutive dimers on the boundary translates 
into a string of n consecutive arrows, all pointing to the left or to the right along the boundary. 
Elsewhere on the boundary, arrows are free to point in any of the three available directions (with 
the only constraint that the full arrow configuration cannot contain loops) . This problem of arrows 



boundary condition on an interval of length 2n, the 



20| ) with arrow configurations, or equivalently spanning 



3 



in the 



with a fixed string of aligned arrows on the boundary has been recently examined in 
context of the Abelian sandpile model. It has been found that the correlation function C2 n behaves 
asymptotically as 

C 2n ^ A n- 1 / 4 e ~ 2nG ^, n large, (5.1) 

where G = 0.915965 is the Catalan constant, and A is a numerical constant. The exponential decay 
is expected and due to the defect of entropy of the n boundary sites which have their arrow frozen. 
On the other hand, the power law n -1 / 4 has been understood within conformal field theory, as the 
correlator of two boundary condition changing fields (and an extra dimension field corresponding 



to the insertion of dissipation, see 



221] ). One of these two fields changes the free arrow boundary 
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condition into say, the right arrow boundary condition, and has dimension —1/8, the other changes 
the right arrow boundary condition back into the free arrow boundary condition, and has dimension 
3/8. The two dimensions add up to 1/4, and account for the exponent in (15. ip . 

Note that these dimensions have been obtained in in the context of the sandpile model 
(equivalently spanning trees), known to correspond to a (logarithmic) conformal field theory with 
central charge c = —2. However, as discussed in the introduction, various results, including those 
of the present article, point to a description of the general monomer-dimer problem in terms of a 
conformal field theory with central charge c = 1, of which the theory with c = —2 would appear 
as a subtheory accounting for the dimer degrees of freedom. In the c = 1 setting, the above power 
law is presumably interpreted as a 2-point function of fields with dimension 1/8. 
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